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PART I. 

INTRODUCTION AND HISTORICAL REVIEW 


The importance of understanding the ionized state of gases 
has become realized in recent years with the enthusiastic attempt to 
solve such problems as nuclear fusion, magnetohydro dynamic generation 
of electricity, propagation of radiowaves in the ionosphere, and ionic 
engine space propulsion. Phenomena associated with gaseous plasma axe 
in general complicated by the presence of a large number of particles 
and the simultaneous occurrence of various events. The study of iso- 
lated basic processes may therefore lead to a basic understanding of 
these phenomena. 

Since the advent. of the microwave method for detecting 

20 

electron densities in ionized gases by S. C. Brown in 1948, the 
study of the basic electron loss processes in afterglows (the time 
following the removal of the ionizing energy source) has received 
considerable attention. Prior to the development of microwave tech- 
niques, electron densities were measured mostly by Langmuir probes, 
which are discussed in detail in various texts on plasma electronics, 

~/r 

e.g., Loeb' 3 . Since the electron collection by the probe disturbs 
the plasma to some extent, it is difficult to interpret probe measure- 
ments. Probe measurements are extremely difficult to make in rapidly 
changing plasma situations, so their main contributions to date are 
results of studies in steady state gas discharges. 

By employing the microwave technique, one is able to measure 
rapid changes in electron density; thus the basic electron loss pro- 
cesses may be studied in the afterglow. The removal of the energy 
source greatly simplifies the experimental conditions. Since highly 
sensitive microwave detection devices are available, signals of ex- 
tremely low amplitude may be used. Measurements may therefore be 
made with probing signals which may be considered to have a negligible 
effect on the phenomena to be studied. 
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Initial studies employing microwave techniques were concen- 
trated on determinations of electron losses by recombination and dif- 
fusion in the inert gases 1 ^ especially helium. It was at first 

assumed that the afterglow contained only atomic ions, electrons, and 

neutral particles, but results of determinations of the recombination 

k 

coefficient in helium were approximately a factor of 10 higher than 

8 

expected. Bates proposed a theory on dissociative recombination of 

the molecular ion, He + , which would explain the higher recombination 

rate, and the presence of the molecular ion in the afterglow was con- 

48 

firmed with a mass spectrometer study by Phelps and Brown . As a 
result of this work, it was realized that afterglow studies by the 
microwave technique must be supplemented by additional experimental 
procedures in order to separate properly the various events which are 
taking place. 

Micro'wave techniques may basically be divided into two groups: 
l) the cavity method, originally introduced by S. C. Brown, in which 
the electron density and collision frequency are related to the shift 
in the resonant frequency of a microwave cavity containing the plasma 
and the absorption loss or "Q" of the cavity and 2 ) the transmission 

p 

method, introduced by Anderson and Goldstein , in which the microwave 
signal is passed through the plasma and the electron density and col- 
lision frequency are related to the phase shift and attenuation of the 
microwave signal. 

Various authors have reported results of measurements in 
nitrogen. In the 1930's, Mitchell and Ridler^ and Tyndall and Pearce‘S 
used a drift tube to measure the mobilities of various gases in nitrogen. 
The first results determined by the microwave technique (cavity) were 
reported by Biondi and Brown^ in 1949* Since that time, other measure- 
ments employing the microwave technique were reported by: Anderson 

and Goldstein^ (1956, transmission); Bialecki and Dougal 1 ^ (1958, 
transmission); Faire and Champion 2 ^ (1959, cavity); Formato and Gilardini 2 ^ 
(i 960 , transmission); and Kasner, Rogers and Biondi ^ ( 1961 , cavity). 

Measurements of electron loss coefficients and cross-sections 
for momentum transfer in nitrogen by the two microwave methods have 
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shown different results. Formato and Gilardini have proposed that 
the measurements by the transmission method might have been carried 
out early in the afterglow before the temperature had relaxed to room 
temperature,, although the measurements were assumed to be carried out 
in a room temperature isothermal gas. Precise measurements of elec- 
tron temperature in the afterglow have not been made to this date, but 
there is no question of the necessity for such measurements in order to 
completely describe the afterglow. 

In this report the nitrogen afterglow has been studied 
using the microwave transmission technique as the basic diagnostic 
tool. A concerted effort has been made to simplify the experimental 
conditions in order to remove some of the questionable assumptions 
made in previous transmission method experiments. A plasma cell 
has been developed which consists essentially of a piece of standard 
waveguide, so that a simple rectangular geometry with well-defined 
boundaries is obtained for the microwave -plasma interaction. A micro- 
wave bridge circuit has been developed so that simultaneous measure- 
ments of the attenuation and phase shift of the microwave signal 
passing through the plasma are obtained. The bridge is calibrated 
directly by a precision attenuator and phase shifter in the microwave 
circuit. 

In addition to the microwave measurements, a newly developed 
gated microwave radiometer was used to measure the electron temperature 
variation in the afterglow. A spectrographic study of the afterglow 
was made and a method for studying the light intensity as a function 
of time is discussed. By this variety of techniques, previous results 
obtained by using the microwave transmission method alone are shown 
to be in error and results of the present experiment are shown to be 
in agreement with cavity measurements. Although some question as to 
the actual recombination mechanism in the nitrogen afterglow still 
exists, the possibilities are discussed and narrowed as much as 
possible. 
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PART II. 
THEORY 


A. Conductivity of a Plasma at Microwave Frequencies 

Considerable work Has been done on the theory of the con- 
ductivity of plasmas at microwave frequencies in the past ten 
years ’ ’ . The permittivity of the medium is usually taken as 

that of free space and the plasma is then completely described by the 
conductivity parameter. Conductivity is defined as the ratio of the 
current density to the electric field intensity producing it, where 
both quantities are assumed to have the time variation where co 

is the angular frequency in radians/second. 

The equation of motion for the electrons is written 1 ^ 


m — + mvV m = qE (l) 

where the second term represents a damping term, and closer inspec- 
tion would reveal that the effective collision frequency for momentum 
transfer -V has been taken independent of velocity. E is the elec- 
tric field of the probing microwave signal and q the electron charge. 
It is advantageous to discuss the derivation of equation 1 to indicate 
the assumptions that have been made. 

The model used is close to that of a free electron gas. It 
is assumed that the heavier ions do not respond to the electric field 
of the probing microwave signal. The gaseous plasma is composed of 
electrons, positive ions and neutral molecules (in the case of an 
electro -negative gas, negative ions would also be present). A dis- 
tribution function would exist for each species of the gas and a 
Boltzmann equation is written for each of the species. The equations 
are normally coupled through the collision terms. Consider the Boltz- 
mann equation for electrons 



c) t 
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where f denotes the distribution function for electrons. The terms 
e 

on the right hand side represent non-explioitly the rate of change of 

the electron distribution function due to collisions where C is the 

ee 

rate of change of the electron distribution function due to collisions 

among electrons, due to collisions between electrons and ions, and 

between electrons and neutral molecules. A great simplification is 

obtained if the gas is weakly ionized so that n g <(<N where n g is the 

electron density and N is the neutral molecule density. It is assumed 

throughout that the plasma is charge neutral, so that the electron 

density is approximately equal to the ion density, n^.. In the weakly 

ionized case, the model assumed is that the electrons find themselves 

in a "sea" of neutrals, so that the electron-electron collisions and 

the electron-ion collisions are negligible in comparison with the 

24 

electron-neutral molecule collisions. Delcroix somewhat arbitrarily 
defines a weakly ionized plasma as one where 


n 


n 


e 

+ N 


<< 10 A 


(3) 


e 

At room temperature and 1 mm of mercury pressure, the mole- 
cular density is approximately 3 x lO 1 ^ particles per cirA Since this 

is in the range of pressures used in the experiments and the maximum 

12 

electron density used in obtaining data was approximately 10 electrons 

3 

per cm , the above inequality holds true. It especially holds true 

in the -late afterglow when electron loss processes are determined, 

9 10 

the electron densities here being in the order of 10 7 - 10 electrons 
per cubic centimeter. 

The equation of motion 1 is obtained by taking the so-called 
second moment of the Boltzmann equation (multiplying each term by the 
momentum mv and integrating over velocity space). The term involving 
the pressure tensor has been omitted, since motion due to pressure 
gradients is negligible with respect to motion caused by the electric 
field in computing the conductivity. The collision term, which is 
normally written c)f/ <9 1 I c, may be expressed as -f 1/ where p 
is the collision frequency for momentum transfer. This assumes that 
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the frequency of the microwave field is sufficiently high so that no 
appreciable change of electron energy or density occurs within a 
single cycle. The distribution function f is then independent of 
time. Strictly speaking, equation 1 should be written 


-Jt (m ^ >} + 


m 


•V <v> = m<a> 
m 


(4) 


where the symbol < > above represents that the quantity has been averaged 
over velocity space. These brackets will be dropped in the following 
discussion. The collision frequency is physically a frequency of 
collision for momentum transfer and has been assumed independent of 
velocity. The acceleration term in equation 4 is a result of the 
electric field intensity, E, of the probing microwave signal, and is 
written as 


a = 


qE 

m 


qE e Jt0t 
o 


m 


(5) 


The solution of equation k becomes 


v = 


qE 

m ( + Jco) 


( 6 ) 


and since the current density j is equal to the charge, q, times the 


particle flow n v 
* e 


n q 2 E 


J = 


m( V + jai) 
m 


(7) 


The conductivity has been defined as the ratio of the current 
density and electric field intensity so that the complex conductivity, 
0“ , may be obtained from equation J. 


<r = 


_v 

m(/V m 


' J®) 


( 8 ) 


This may be separated into real and imaginary parts, and <77, 
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where 


2 

_ _ n e* 

f .2 2s 

ml <y + co ) 
' m 


( 9 ) 


and 


<r. = 

i 


2 

-n q co 


m( 2 + co 2 ) 
v m ' 


( 10 ) 


It is interesting to note here that in the absence of 
collisions ( -V = 0), the conductivity is purely imaginary and the 
model of the free electron gas is obtained. 

The assumptions made in deriving the conductivity are valid 
except for the questionable assumption that the collision frequency 
is independent of velocity. The form indicated in equation 8 is in 
wide use, but in situations where collisions are important, e.g., 
in obtaining collision cross-sections, equation 8 must be taken as a 
first approximation. As it turns out, in the late afterglow at low 
pressures, the effect, of collisions is negligible and electron loss 
process coefficients obtained after using equation 8 are valid. 


B. Electron Losses in a Decaying Plasma 


1. Introduction 

The electrons formed in a gas discharge will disappear by 
three basic processes, when the source is removed. First, the elec- 
trons will recombine with positive ions. This process will be called 
volume recombination, or simply recombination, in the ensuing dis- 
cussion. Secondly, the electrons will diffuse to the wall and they 
will disappear by "wall recombination". The third process is only 
important in the case of electro-negative gases, and that is the 
attachment of electrons to neutral molecules forming negative ions. 

The type of diffusion encountered is ambipolar. The meaning 
of this term will be explained in the next section. The loss mechanisms 
are first discussed singly as though only one mechanism were taking 


7 



place at a time. After that the more difficult problem of combinations 
of loss mechanisms is discussed. 

2. Arobipolar Diffusion 

a. Positive ion-electron . Assume that the mean free path 
is much smaller than the dimensions of the container. Then normal 
diffusion of particles occurs as a result of a density gradient and 
is characterized by 

P = -D Vn (ll) 

where P is the particle flux or flow and is equal to the product 
nv ; and D is the diffusion coefficient and has the dimensions of area/ 
time. Electrons and ions recombine readily at the walls of the con- 
tainer, since the molecules or atoms of the wall material are available 
as third bodies to take up the liberated energy. The density of 
charged particles is thus essentially zero at the walls. This results 
in a density gradient and electrons and positive ions will diffuse to 
the wall. Because of the difference in masses, the electrons would 
diffuse faster than the positive ions. This phenomenon leads to a 
charge separation resulting in a space charge field, which is written 

as E . We must therefore conclude that the diffusion of the elec-' 

S 

trons is retarded by the space charge field, so that both types of 

particles diffuse at the same rate. The simultaneous diffusion of 

both species may be characterized by an ambipolar diffusion coefficient, 

D , which may be found as a function of the diffusion coefficients 
8 # 

and mobilities of the electrons and positive ions. On the assumption 
that the charge carriers consist of electrons and one type of posi- 
tive ion, modified flow equations (equation 11 ) may be written as 
follows: 


P = -D \7n - n u E 

e e v e e'e s 

r + = -D + V n + + n + h + E s 


( 12 ) 


8 



where jj. is the mobility (velocity per unit field). The second term 

on the right hand side describes the flow due to the space charge 

field E . The densities of the electrons and ions are n and n 
s e + 

respectively. The basic assumption of a charge neutral plasma 
(n g = n + )and the assumption of equal flow rates lead to the pro- 
perty: 

f - r = r a 3 ) 

*r e 


The space charge field and P may be found from equations 12 with 
the result 


E 


s 



+ 


V n 

n 


( 14 ) 


where Vn/n may refer to either electron densities or ion densities 
and 


p < * ^+ E e ) 

( + •*+ ) 

Comparison of equations'll and 15 
coefficient is 


V n (15) 

shows that the ambipolar diffusion 


D 

a 


^ e D + + ^ + D e 


(16) 


This expression may be simplified if the distribution func- 
tion of both the electrons and positive ions is Maxwellian. In such a 

situation, the mobilities and diffusion coefficients are related by 

Ij. 

the Einstein relations 


D. 


kT 


u 

J_e 

D 


kT 


(17) 

rv-j. u jxx 

+ g e e 

where k is the Boltzmann constant and q is the absolute value of the 
charge on an electron. The ion temperature will be essentially equal 
to the neutral gas temperature and has been written as T . With the 
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knowledge that the electron mobility is much larger than the ion mo- 
bility, the ambipolar diffusion coefficient becomes 



and in the case of an isothermal plasma where the electron temperature 


is also equal to the neutral gas temperature 


kT 

D 2D = 2|i, ( 1 - 9 ) 

It is seen, then, that if certain conditions are fulfilled, obtaining 

the ambipolar diffusion coefficient leads to knowledge of the posi- 

1 

tive ion diffusion coefficient and the mobility of the ion. 

b. Electro-negative gases . Since gases composing the at- 
mosphere are of special interest, it would be assumed that considerable 
work would be published on properties of oxygen. Neutral oxygen mole- 
cules have an affinity for electrons, thus forming negative ions. It 
is difficult to interpret measurements in such an electro -negative gas; 
thus, only preliminary results have been published to date. 

For the study of ambipolar diffusion in an electro -negative 
gas, assume that the charge carriers consist of one type of positive 
ion of density n , one type of negative ion of density n_, and elec- 
trons with density n . 

The flow equations in this case become 



-D 

e 





7n, + nu E 
v + +^+ s 


r = -D yn - npE s 

The neutrality condition of the plasma requires that 


and 


n = n + n 
+ e 



( 20 ) 


( 21 ) 
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It is assumed that the diffusion has been in effect long 
enough so that the spatial distribution of carriers is characterized 
by the fundamental diffusion mode. More will be said about this later, 
but in effect the assumption required is that the spatial distribution 
satisfying the breakdown mechanism has disappeared and has been re- 
placed by one which satisfies the diffusion of the particles. If this 
condition is satisfied, 


F n e V n 

n n n 

e + 


( 22 ) 


Under these conditions, equations 20 may be solved for the space 
charge field, 


E 

s 


Dn + D n - Dn 

e e - - + + 

n, u + n ll + n la 

+ i~ + g r- e _r-_ 



(23) 


in which the last factor in n may be replaced by any one of the three 
factors in equation 22. 

If the flow equations are now written in terms of an ambi- 
polar diffusion coefficient which accounts for the space charge effect: 


-D Vn 

ae e 

-D a+ Vn + (210 

-V P n - 

The ambipolar diffusion coefficients may be found to be 



r = 

e 

P + = 

r = 


11 



where a 


n 


n 

e 


or 


(l + a) = 


n 


n 


In this case the ambipolar diffusion coefficients are 
different for each type of charge carrier and are dependent on the 
relative concentrations of the various densities. Note that if the 
negative ion density disappears (n = a = 0), then D = D and 

“ 8fr 8,0 

these quantities are equal to the amhipolar diffusion coefficient 
found in the previous section. 

These formidable expressions may be simplified, provided 
that the negative ion density is not large compared to the electron 
density. The mobility of the electron is larger than the mobilities 
of the massive positive and negative ions, so that the denominators 
of equation 25 may be written as approximately |i g , provided that "a", 
which is the ratio n_/n e , is not too large. It is not difficult to 
realize this situtation in practice. Under this condition and using 
the Einstein relations (equation 17 ), the expressions 25 become 



For an isothermal plasma (T g = T ), these expressions reduce 


further to 


D = 2 (l + a) D 
ae v ' + 


D = 2 (1 + a) 
8,~ 


D + D - 

D 


( 27 ) 


D = 2 
a+ + 


It is seen in equations 27 that the ambipolar diffusion 
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coefficient of the positive ions is little affected by the presence 

of the negative ions in an isothermal plasma. The ambipolar diffusion 

of the electrons is increased by the ratio of the densities of the 

positive ions and electrons. D is smaller than either D or D , 

since the ratio D /D is less than unity. In the absence of nega- 
- e 

tive ions (a = 0), the result of the previous section is obtained 

with D = D = 2D . It must be remembered that the assumption 
ae a+ + 

has been made that the number density of the negative ions is not 
too large with respect to the density of electrons. In a situtation 
where this is not true, the more general expressions 25 would have to 
be used. 

3* Electron Loss by Diffusion 

The flow equation 11 is written with a free diffusion 

coefficient or ambipolar diffusion coefficient, depending on the par- 

h 

ticular experimental situation. Allis has discussed the requirements 
in some detail. As might be expected, the diffusion will be ambi- 
polar, provided the characteristic Debye shielding length is less than 
the dimensions of the vessel. The Debye length (rationalized MKS units 
are used) is given by 

(28) 

» n q 
e 

Q 

At room temperature and at an electron density of 10 per 
cubic centimeter, which is lower than the minimum detectable density 
in the experiments, the Debye length is less than a millimeter. This 
is much smaller than the dimensions of the vessel, so ambipolar dif- 
fusion is expected. 

The continuity equation, which is obtainable from the first 
moment of the Boltzmann equation, is 

+ 7. P = 0 (29) 

Inserting the equation of flow in terms of the ambipolar" diffusion 
coefficient leads to the equation describing the time and spatial 
variation of the density. 
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( 30 ) 


J n 
c) t 


= D V 2 n 

a r 


It is assumed here that one type of positive ion and electrons are 
present and that n = n = n. D then is given by equation 1 6. 

0i £1 

Equation 30 is easily solved, when the density is subject to simple 
boundary conditions, by a separation of variables technique. Two 
specific problems of interest will be worked out-- first, diffusion 
in a rectangular box and secondly, diffusion in a right circular 
cylinder. 

The plasma cell in the experimental arrangement consists 
of a length of rectangular X-band waveguide. With a coordinate system 
in such a geometry as indicated in Figure 1, the solution to equation 
30 is 


oO 


n(x,y,z,t) = 

k,J?,m=l 


m 


exp 


D t 
a 


a 


|sin(^)sin(^)sin(S^p) (31) 


k,J?,mJ 


(|i) 2 + ( Mf + (2*) 2 

j\_ 0 is known as the characteristic diffusion length, 

a, b, and c are the dimensions of the rectangular vessel. The boundary 
conditions are that the electron density is zero at the walls because 
of wall recombination. 

The coefficients K, . are the Fourier coefficients; they 

, m 

depend on the initial electron density distribution and are given by 
c b a 

n(x,y, z, 0 )sin(^~ )sin(4^)sin(2^-) dx dy dz ( 32 ) 

000 


K 


'k,i,m 



where 


A 


k, 9,m 


A 
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Fig. 1. Rectangular geometry for diffusion problem. 



Fig. 2. Cylindrical geometry for diffusion problem. 



It is noticed in equation 31 that if one waits long enough, 

the fundamental diffusion mode (with k = J( = m = l) will predominate 

and a semi logarithmic plot of n vs. t will he a straight line. The 

slope of the plot will give the amhipolar diffusion coefficient, D , 

s> 

since -A . , .depends on the dimensions in a known fashion. 

Another specific problem of interest is diffusion in a 
right circular cylinder. This problem has application in glass tube 
vessels and in cylindrical cavities. Assume the coordinate system 
origin to be at the bottom center of the cylinder as shown in Figure 
2 , and axial symmetry for the density function. The solution of 
equation 30 is then 


n(r,z,t) = 


i, 0 v 0 


( 33 ) 


where 


A 


= (^) 2 + (£l ) 2 
M 


is the kth root of the zero order Bessel function 1 . 
k 0 

The Fourier -Bessel coefficient, K^. ^ , may be found from the initial 
electron density distribution n(r, z,o) and is given by 



R L 


o o 


n(r,z,o)rJ Q (A k |)sin(| iI z) 


dz dr 


( 34 ) 


As in the case of rectangular boxes, the ambipolar diffusion 
coefficient may be found from the slope of the semi logarithmic plot 
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of the electron density vs. time when the fundamental diffusion mode 
is predominant. 


k . Electron Loss by Recombination 

Recombination between positive ions and electrons may take 

place in a variety of ways, i.e., radiative, di -electronic, dissociative, 

bb 

and three -body recombination. Massey and Burhop discuss the different 
processes in some detail, so they will not be recounted here. Each 
one of the processes may have associated with it a coefficient of 
recombination which depends on pressure and . temperature in a particular 
way. The microwave method, measuring the electron density as a function 
of time, leads to a recombination coefficient which is related to more 
than one recombination process. In many physical situations, however, 
the recombination coefficient for one process predominates and may 
be recognizable from its pressure and/or temperature dependence. 

Assume that the plasma exhibits the property of quasi - 
neutrality (n + = n = n). Then, because of recombination, 


dn 

dt 


dn. 


dn 


dt 


dt 


e 2 

= -an , n = -an 
+ e 


(35) 


The quantity a is the recombination coefficient, 
solution of equation 35 is 


The 


n(t) 


n 


Tzy 


+ at 


(36) 


In the absence of other loss mechanisms, the recombination 
coefficient may be obtained from the slope of a linear plot of l/n 
as a function of time. 

The relative loss rate (l/n)dn/dt is proportional to the 
electron density. Volume recombination (remote from the walls of 
the container) would make the density uniform, since the relative 
loss rate is greater in regions of higher concentrations. 

As mentioned in part B 2 a, electrons and ions recombine 
readily at the wall of the containing vessel, so that the concentration 
at the wall may be assumed zero. Therefore, there is always a loss by 
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diffusion to the walls, hut situations exist where this loss may he 

negligible when compared to volume recombination. Since the rate of 

2 

electron loss by volume recombination is proportional to n , recom- 
bination would predominate with a high degree of ionization or at 

high pressures where the diffusion loss is reduced. 

12 

Gray and Kerr have done extensive work examining the re- 
liability of results for the recombination coefficient based on l/n 
vs. time plots. They have made computations which include the loss by 
diffusion, and have concluded that a linear plot of l/n vs. t does not 
necessarily imply a loss by recombination, unless the range of electron 
densities covers at least an order of magnitude. Some of the earlier 
recombination coefficients reported by investigators who used the 
microwave method may be suspect on the basis of Gray and Kerr's work, 
since linear plots of l/n vs. t were employed in which the electron 
density ranged over a factor of only two or three. In fact, Kerr, 
employing the equation describing loss by diffusion alone, has pointed 
out that plots of l/n vs. t may appear to be linear over small ranges 
of electron density. Extreme caution is advised when such plots are 
to be interpreted. 


5. Electron Loss by Attachment 

Electrons may be lost by attachment to a neutral atom or 

molecule whose outer electronic shell is nearly filled. The most 

common example is oxygen, which shows a high affinity for electrons. 

In the early days, investigators measured electron attachment rates 

by allowing electrons to drift through an attaching gas, separating 

electron and negative ion currents by various means and comparing the 

6 22 17 

magnitudes of the electron and negative ion current ' ’ J . 

When attachment is the predominant loss of electrons, the 

microwave method can be used to measure the attachment coefficient, 

a) . The rate of electron loss may be written as 
r a 



(37) 


where n is the electron density. 
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The solution of equation 37 is simply 

- Vt 

n(t ) = n(o ) e a (38) 

The attachment coefficient can he obtained from the slope of a semi- 
logarithmic plot of electron density as a function of time; the slope 
should increase linearly with pressure at constant temperature, since 
the rate of attachment is proportional to the number of neutral par- 
ticles present. 

The relative loss rate, ( /n) dn /dt, is equal to a constant 
so that the spatial distribution of electrons is unaffected by the 
attachment process. 

6. Electron Loss by a Combination of Mechanisms 

The electron loss processes in a plasma do not occur singly, 
in general, but processes occur simultaneously. As we shall see 
shortly, this leads to non-linear differential equations and it is 
difficult to obtain the individual loss coefficients. In some physical 
situations, one process may predominate so that linear approximations 
may be made. In these situations the individual loss coefficients 
may be obtained. 

Two examples will be discussed. The first case is in a gas 
in which attachment does not occur. The remaining two loss mechanisms., 
volume recombination and ambipolar diffusion, occur simultaneously. 

In the second case, that of attaching gases, the equations are written 
for simultaneous loss by volume recombination, ambipolar diffusion 
and attachment. Procedures are indicated which allow solution under 
specified physical conditions. 

a. Loss by recombination and ambipolar diffusion . Assume 
that one type of positive ion and electrons are present in the after- 
glow, and that the positive ion density is equal to the electron density. 
The equation for the rate of electron density loss is 

= D & v 2n " (39). 
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where n is the electron or positive ion density, D is the ambi polar 

El 

diffusion coefficient given in section 2 a, and a, is the recombination 
coefficient. The boundary condition is that n is zero at the walls of 
the container, because of wall recombination. The equation is non- 
linear, separation of variables technique fails, and no exact solution 
is known. 

The key to the method of attack in this problem is the re- 

2 

combination term Chi . Since the loss by recombination is proportional 
to the square of the electron density, this term diminishes rapidly so 
that in the late afterglow the principal loss is by diffusion. One 
would expect, then, that the spatial distribution would approach, that 
of the fundamental diffusion mode as described in part II 3. For 
example, in the case of rectangular vessels as indicated in Figure 1, 
the spatial distribution will become sinusoidal, as shown by equation 
31. Once the spatial distribution is known, equation 39 is reducible 
to an equation in one variable, time. 

The effect of the recombination term and higher order dif- 
fusion modes can be shown more rigorously by a machine solution of 

47 

the finite difference form of equation 39- Oskam has done this 
assuming diffusion in one dimension, using a uniform spatial distri- 
bution as the initial condition. His results indicate, as expected, 
that the. loss by higher order diffusion modes is negligible after a 
sufficient time has elapsed in the afterglow. After this critical 
time, the loss can be described essentially by recombination and the 
fundamental diffusion mode whose spatial density distribution is 
known. This critical time occurs when the density is a factor of 
four or five above the value obtained from the straight line repre- 
senting the decay of the fundamental diffusion mode alone. 

If the sinusoidal spatial density distribution of the funda- 

•X* 

mental diffusion mode is assumed, then equation 39 reduces to 
dn (t) 

~dt = “ kn o (t) ' ^o (t) (40) 

where n Q is the electron density at the center of the rectangular 

* See Appendix 
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vessel whose dimensions are a x To x c, and the constant is given "by 


, _ -p, 2 / 1 , 1 , 1 \ 

k “ D a n ^ 2 , 2 2 


(41) 


The solution of equation 40 is 

-kt 


n 0 (t) 


C^ke 


1 - aC-^e 


-kt 


n (o) 
o 

where C, = — — — — 7 — r- 

1 k + an (o) 


(42) 


Obtaining the ambipolar diffusion coefficient and recombi- 
nation coefficient from a plot of n (t) as a function of time and 
equation 42 is somewhat of a problem in itself. 

Some authors 1 ’ write . equation 42 in the form 


1 + 


n Q (t) 


Qh o (t) 


1 + 


r^o) 

an (o) 

o ' 


k 


-kt 

e 


(43) 


and adjust the factors a, k, and n (o) for a best fit to the data 
points (comparing both sides of equation 43). In the present work 
a different approach is used. Equation 40 is written in the finite 
difference form 



A n Q (t) 

A t 


k - an (t) 
0 


(44) 


and the left hand side of the equation is plotted as a function of 

n (t). The coefficients k and a are found from the intercept and 
o 

slope of the resulting straight line, using the method of least 
square deviation. In order to achieve accuracy with this method, n Q (t) 
must be found at small increments of time, which requires a large 
amount of data and routine calculation. Since the routine calculation 
could be done by machine (IBM 6 50 Digital Computer was used), it was 
felt that this method was superior to the method employing equation 
43. The present method leads directly to D and a, while the previous 
method required adjustment of parameters for a best fit to data points. 


21 



b. Loss by recombination, ambipolar diffusion, and attach- 
ment . The equations describing the charged particles of the plasma 
in the afterglow of an electronegative gas are: 


^ n e 

d n _ 


p 

= D V n - -17 n 

ae ' e v b. e 

p 

= D 17 n + ^ n 

a- r a e 

2 

= D , V n , - CC n, n 
a+ Y + - + - 


a n n, 
e e + 


a n + n 


a n, n 
e + e 


(^ 5 ) 


The electron density, positive ion density, and negative 
ion density are n , n , and n respectively. The first term on the 
right hand side of each equation represents the loss by ambipolar 
diffusion and the ambipolar diffusion coefficients are given by 
equations 25 • Af is the attachment coefficient, a is the negative 
ion-positive ion recombination coefficient, and cn e is the electron- 
positive ion recombination coefficient. The first equation relates 
the rate of change of electron density to losses by ambipolar dif- 
fusion, attachment to neutral particles, and recombination with posi- 
tive ions. The second equation relates the rate of change of the 
negative ion density to losses by ambipolar diffusion and recombi- 
nation with positive ions. This equation has a source term 1 / n 

8 . © 

which describes the formation of negative ions by attachment. The 
last equation relates the rate of change of the positive ion density 
to losses by ambipolar diffusion, and recombination with both electrons 
and negative ions. 

Since the microwave method of measurement is responsive 
only to the electron constituent of a plasma, only the electron den- 
sity as a function of time may be obtained. The possibility of ob- 
taining electron loss coefficients from equations 45 by the microwave 
method therefore appears remote, expecially when it is noticed that 
the equations are not only coupled through the various constituent 
densities, but also through the ambipolar diffusion coefficients (in 
equations 25> the factor "a" is the ratio of n to n^). There axe 
some simplifications which may be made in certain physical situations, 
however . 
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The easiest interpretation of data occurs when the dif- 
fusion terms in equations 45 are negligible. Such a situation occurs 
at higher pressures, where the diffusion loss is hampered by the in- 
creased elastic collision process with neutral particles. If, in 
addition, the negative-ion constituent of the plasma is negligible, 
the requirement of charge neutrality requires that n g = n + and the 
equation describing the electrons in 45 becomes uncoupled from the 
remaining two. 


^ n e 


yjj n 
a e 


a n 
e e 


(1*6) 


This equation may be used to obtain the attachment coefficient and 
recombination coefficient by a method similar to the one described 
in the previous section. 

The negative ion density may become negligible for two 
different reasons. First, the cross-section for recombination of 
the negative ions and positive ions is much larger than the cross- 
section for the recombination of the electrons and positive ions. 

The negative ions will therefore disappear much faster than the 
electrons. Secondly, the production of the negative ions may be re- 
duced by using low partial pressures of the electronegative gas and 
using a "buffer gas" to increase the total pressure to a value where 

34 

loss by diffusion is negligible. Kasner, Rogers and Biondi have 
used this method in the study of recombination in oxygen. Neon was 
used for a buffer gas. The ionization potential of neon is suffi- 
ciently high so that one observes only ions of oxygen. Electron 
attachment should not be a significant factor since the oxygen 
partial pressure was low compared to the total pressure. In this 
last case, equation 46 would be reduced to loss by recombination only 
as discussed in part II B 4. 

In the presence of ambipolar diffusion the situation appears 
to be extremely difficult. Equations 45 then are nonlinear and no 
general solutions are known. The densities are functions of both 
space and time and the spatial distributions are difficult to ascer- 
tain. At low electron densities, with a low degree of ionization, 
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the electron-positive ion recombination loss will be small, and. the 

spatial distribution will approach that of the lowest diffusion mode 

(described in part II B 3)* Even though the spatial distribution is 

known, and recombination is negligible, equations 45 are still coupled 

and additional assumptions on the presence of the negative ion are 

required before the equations are soluble. If the negative ions 

formed disappear very rapidly, D becomes independent of the den- 

8,6 

sities, so that the first of the equations 45 is decoupled from the 
other two and is soluble. 

The problem in the general solution. p>f equations 45 is 
quite insurmountable, and, as pointed out in the previous discussion, 
only in a few special cases may the loss process coefficients be ob- 
tained by microwave methods in electronegative gases. In these special 
cases, a large amount of ingenuity is required on the part of .the experi- 
menter in setting up physical conditions so that the . loss process 
coefficients may be interpreted from the data. 

C. Relation of Microwave Measurements to Plasma Parameters 

The complex conductivity of a weakly ionized plasma at 
microwave frequencies has been found in part II A. If the complex 
conductivity is obtained experimentally, then the electron density and 
collision frequency for momentum transfer may be found by equations 
9 and 10. It is now necessary to relate the measurable quantities 
of phase shift and attenuation of the microwave signal to the complex 
conductivity. Two cases are of interest. If recombination is pre- 
dominant in an afterglow, the spatial distribution of the electrons 
would be expected to be uniform (see section on recombination, II 
B 4). In this case, the relation of phase shift and attenuation of 
a microwave signal passing through the plasma to the electron density 
and collision frequency is easily derived. 

If diffusion is the predominant loss mechanism, the electron 

gas is not uniformly distributed and the relationship of the measurable 

quantities and the plasma parameters is not straightforward. In this 

second case, an excellent approximation method has been developed by 

2 

Anderson and Goldstein . This method is briefly outlined, and the 
simplifying assumptions required for solution are identified. 
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1. Uniform Electron Density Distribution 

The wave equation for the electric field intensity, E, 
of an electromagnetic wave passing through a homogeneous lossy medium 
is well known and is given hy 

y 2 il = joop ( a" + jtoe) E (47) 


It has heen assumed here that the time variation of E is 
sinusoidal and that the time variation is written in the form e^ 03 ^. 
cd is the radian frequency and the medium is described by p, e, and 
which are the permeability, permittivity, and conductivity res- 
pectively. The permeability and permittivity of free space will be 
used, and the plasma will be described by a complex conductivity 
as discussed in part II A. 

The solutions of the wave equation may be divided into the 
familiar transverse electric and transverse magnetic modes in a wave- 
guide. The following discussion assumes the dominant TE 1Q mode, 
since the experiments were carried out in a rectangular waveguide 
at a frequency where all other modes were cut off. A parallel dis- 
cussion may be used for a cylindrical waveguide. If the geometry 
is assumed to be that of Figure 1, with field variations in the z 
direction of the form e , the complex propagation constant,"^, 
may be found by solving the wave equation with the boundary conditions 
of Figure 1. For the TE^q mode, Y is 


Y = + ^o^ = a+ ( 48 ) 

The real part of the propagation constant, a, is the 
attenuation/length, and the imaginary part, P, is the phase shift/ 
length. For a uniform plasma filling the cell, the conductivity is 
given by equation 8, which is rewritten in the form 




^m - J 

V 2 + CD 2 

m J 


(49) 
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It is noted from equations 48 and 49 that if the complex 
propagation constant is measured, then the electron density, n, and 
collision frequency, , are obtainable. 

If Z is the length of the plasma filled waveguide cell, 

-qcJL 

an electromagnetic wave is attenuated by a factor e in passing 
through the plasma. (3 j2 is the phase change of the wave as it passes 

through the cell. The phase shift actually measured is the change in 
phase that takes place when the cell is changed from the un-ionized 
condition to an ionized condition. This may be written as 

<#>=e c i -fU (50) 

where (f) is the measured change in phase, is the phase shift 
constant when there is no plasma in the cell ( CT~ = 0) and (3 is the 
phase shift constant when the gas is ionized (given by equation 48). 

The method of experimentally determining <j> and a Jl will be des- 
cribed in part III in the description of the microwave bridge circuit. 

Equation 48 may be separated into real and imaginary parts, 
so that explicit equations may be found for the electron density and 
collision frequency in terms of the phase shift and attenuation. The 
results are given conveniently in equation 51 in terms of dimension- 
less variables which are then defined. 
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where 
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0 ) is the plasm^. frequency, cd is the operating radian frequency, 

P — l/? 

c is the velocity of light (m- 0 € q ) > an & (x> c is cut-off radian 

frequency for the TE^q mode for the empty waveguide (equal to jtc/a). 

2. An Approximate Method for Non-Uniform Electron Density Distributions 

It was determined in part II B 4 that in an afterglow con- 
trolled by recombination, the electron density approaches a uniform 
spatial distribution. In such a situation the attenuation and phase 
constants of an electromagnetic wave passing through the plasma may 
be calculated exactly by the method described in the previous section. 
When the plasma is non-uniform, it is difficult to relate the measurable 
phase shift and attenuation to the electron density and collision fre- 
quency of a plasma. In a waveguide containing a homogeneous medium, 
the solution of the wave equation separates into modes which are inde- 
pendent of each other. In the inhomogeneous case, it is very diffi- 
cult to obtain an exact solution that may be used in practical measure- 
ments. An excellent discussion of an exact method of solution for 
waves on inhomogeneous cylindrical structures has been given by Adler 1 . 

The method used here is an approximate one, but one which 
readily adapts itself to the calculation of electron density and col- 
lision frequency of a plasma from experimentally determined phase shift 
and attenuation measurements. The basic assumption is that the funda- 
mental mode propagated in the empty waveguide is perturbed only by a 
small amount when a weakly ionized plasma fills the waveguide. This 
assumption is valid at high frequencies, where the conductivity, 0~ , 
in equation 48 is small. This method is not applicable in situations 
where the plasma frequency approaches the value of the operating 
frequency, where the conductivity term becomes appreciable. 

Maxwell's equations for the plasma medium may be written 
in terms of the complex conductivity as 


V x E = - J (0 n Q H 


v x H = ( 0“ + j 0 ~. + j CD e ) E 

y r 1 o 


(52) 
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If E and H are expressed as complex values, the time averaged complex 
power flow through a closed surface is given hy 


Power = 


X 


l/2 (E x H*) • da 


( 53 ) 


where the asterisk denotes the complex conjugate. By use of the 
divergence theorem and Maxwell's equations, this may he manipulated 
into the form 


X 


(E x H ) 


da = 




- — -sf \ . 0 i * - 

CTE*E dv + jco \ (e + ) E*E dv 

r 0 0 U) 

v 


joD(j. Q H-H* dv 

v 


(54) 


Only the ends of the plasma filled Waveguide section con- 
tribute to the integral on the left hand side of equation 54, since 
the walls are assumed perfeatly conducting. The field intensities 
may he expressed as the vector sum of two quantities; the first a 

component in the direction of propagation (H and E ), and the second 

z z 

a vector transverse to the direction of propagation (H and E. ) . For 
TE modes the field intensities are related hy the expressions 


x E , = 
z t 


j^M. c 

T“ 


Hi 


and 


H t = 


- r 

T 2 + 


2 

co u e 
*0 o 


Vz 


where y H is the gradient in the transverse plane. Thus all field 
"0 z 

quantities may he expressed in terms of the transverse electric field 
intensity, E . If we assume field variations in the z direction along 

"t ** y 

the waveguide axis of the form e * where yf = a + jp, then we may 
divide equation 54 into real and imaginary parts with the result: 
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where 

= phase shift constant (radians/meter) of a free space 
(unguided) wave; it is equal to co/c. 

j 3 = phase shift constant of the guided wave in a plasma- 
filled waveguide. 

= phase shift constant of the guided wave in an empty 
waveguide. 

The integrations are to he carried out over the cross-section 
of the waveguide. If the electrons are distributed uniformly, the 
conductivity terms can he taken outside the integral signs and the 
integrals would cancel. This leads to the result given in the previous 
section. Equations 55 may he used to integrate the non-uniformities 
of the conductivity over the cross-section of the waveguide and it is 
noticed that, in effect, the conductivity is weighted and averaged hy 
the square of the electric field intensity. 

If is the length of the plasma-filled section, equations 
55 may he written in terms of the total attenuation a $■ and the phase 
shift ^ which are measurable quantities. <p again is the change in 
phase shift that takes place when the cell is changed from the un-ionized 
condition to the ionized condition and is given hy equation 50- 
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The special example of interest here is to relate the elec- 
tron density and collision frequency to the measurable attenuation and 
phase shift when diffusion is the predominant loss mechanism. It was 
shown in part II B 3 that in the rectangular waveguide cell, diffusion 
leads to a sinusoidal spatial distribution of electrons. The collision 
frequency for momentum transfer, will be assumed to be spatially 

independent. This should be approximately correct, since is a 
function of the electron energy, and it is assumed that thermal gradients 
are absent in the late afterglow. Under these assumptions, . ( 5 “ and 
<TT both have sinusoidal variations in both the x and y directions. 

E for a TE mode has only a y direction component and the amplitude 
0 -L*. 

varies sinusoidally in x. Using these variations and letting n Q be the 
electron density along the centerline of the waveguide cell, equations 
5 6 become 



Finally.} equations 55 and 57 may be solved for the electron 
density and collision frequency. The results, written in dimension- 


less form, are 




30 



where the parameters are those defined in equations 51. The same 
phase shift and attenuation measurement would lead to an electron 
density at the center of the waveguide for the sinusoidal distribution 
3it 2 /l6 (approximately 1*85) times the value obtained in a uniform 
distribution. 

Equations 51 are used to obtain the plasma parameters in 
the early afterglow when the electron density is high and recombi- 
nation is the principal loss mechanism. Equations 58 are used in the 
late afterglow when the plasma is controlled by diffusion and the 
spatial distribution is sinusoidal. 

D. Measurements of the Radiation Temperature 

The problem of electron temperature measurements during the 
plasma afterglow has to date received little attention. Originally 
investigators estimated the time required for electron energy to 
relax to thermal equilibrium with the heavier gas constituents fol- 
lowing the removal of the excitation source. They assumed that the 
time required to achieve thermal equilibrium is dependent primarily 
on the electron energy decrease through elastic collisions of the 
electrons with neutral particles. Results of different microwave 
measurements in nitrogen afterglows have been in disagreement. 

Formato and Gilardini 2 ^ have pointed out the possibility of the pre- 
sence of energetic electrons coming from ionizing collisions of long- 
lived metastable molecules, providing for an increase in electron 
energy. Thus time estimates for reaching thermal equilibrium based 
on elastic collisions may be too short, and some results of microwave 
measurements reported as carried out with room temperature electrons 
may have, actually involved more energetic electrons. Thus the need 
for electron temperature monitoring is established as one of the 
fundamental requirements for a complete interpretation of afterglow 
experiments . 

The basic problem of relating the electron temperature to 
the power radiated at microwave frequencies is in accounting for the 
absorption. In a model of a radiating hot body, Kirchhoff 's radiation 
law relates the temperature of the body to the ratio of emission to 
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absorption of the body. Bekefi 10 '’ 11 has used this classical model 
to describe a method of measuring electron temperatures without know- 
ledge of the absorption of a plasma. Kirchhoff 's law may be written 
as 



(o>,T) 


where P^dco is the radiation intensity in the frequency range cd to 
cd + dco, B (<jd,T) is the intensity of black-body radiation, and is 
the absorptivity of the plasma. A varies from unity for a black body 
to zero for a non-absorbing body. 

In order to determine the temperature of a body by direct 
measurement of the radiation intensity, P , the absorptivity, A , 
would have to be obtained separately. The technique used in this 
thesis is to compare the radiation from the plasma with a known variable 
standard radiating through the plasma to the observer, as shown in 
Figure 3* 


Let the plasma temperature be denoted by T and the tem- 

P 

perature of a known variable black -body source as T . The emission 
from the known source is B (co,T ). If this radiation passes through 

3 

the plasma, a fraction A^ is absorbed. Referring to Figure 3> an 
observer on the opposite side of the plasma would therefore receive 
a fraction 1 - A of the standard source’s radiation intensity. Adding 
the intensity of the plasma radiation, the total intensity seen by an 
observer on the opposite side of the plasma from the variable standard 
would be 


P A , = (l - A ) B (cd,T ) + A B (cd,T ) (60) 

01 v CD 7 ' 7 S 7 CD v 7 p 7 

The reflectivity is assumed to be zero. If the plasma 
were removed from the path between the observer and the known source , 
the total intensity received would be simply 


P Q2 = B (cd,T s ) (6l) 

The difference in intensities in the two cases is 
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Arrangement for radiation measurements 




When ,/IFq is zero, = T g , regardless of the value of the 
absorption term A • 

The measurement of the plasma temperature is obtained by 
periodically inserting the plasma between the observer and the stan- 
dard, adjusting the standard until ^ Pq = 0. When this null is 
obtained, the plasma temperature is equal to the known standard tem- 
perature. The details of the experimental arrangement for temperature 
measurements and a description of the radiometer are included in 
part III D. 

Since the emission at microwave frequencies in weakly ionized 

plasmas is mainly the result of the motion of free electrons (brems- 

strahlung emission from collisions of electrons and neutrals), the 

temperature measured will be that of the electron gas. However, in 

assigning an electron temperature, the velocity distribution is assumed 

Maxwellian. If the velocity distribution function is other than 

Maxwellian, the parameter measured is a fictitious temperature, still 

related to the mean electron energy, but through the velocity d±s- 

12 

tribution function, which may not be known . A Maxwellian distri- 
bution is expected in the afterglow unless a source mechanism, such 
as the release of metastable energy, is present. 


E. Relation of Cross-Section to Collision Frequency 
Collision frequency in a gas is defined as the mean velocity 
of a particle divided by the mean distance traveled between collisions. 
In a general sense, the word collision may describe any interaction 
such as recombination, excitation, "billiard ball" type elastic col- 
lision, etc. In a weakly ionized gas the dominant collision process 
involves the elastic scatter of electrons by neutral particles. It 
is advantageous to introduce a cross-section for momentum transfer 
to describe such a process. 

40 

The total scatter cross-section is defined as: 
n 
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(63) 



where I ( 0 ) is the angular distribution function for scatter into 
the angle between Q and Q + d£ , where 0 is measured in center 
of mass coordinates. 

For an elastic collision between an electron and a neutral 

lj-1 

molecule, the speed of the electron is essentially unchanged . The 
change in momentum of the electron along its initial direction is thus 

Change in momentum = |j,v (l - cos0) (64) 


where (i is the reduced mass 


mM 

^ m + M 


(65) 


and the neutral particle is considered to be stationary. 

The mean fractional change in momentum is therefore related 
to the mean value of (l - cos 6> ), which is 


^1 - COS = 
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where Q is 
mi 


Tn 
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- 2n (l - cosd ) I ( & ) sin© d© 

O 0 


(67) 


and Q is referred to as the cross-section for momentum transfer. 


m 


If the model of the collision is taken as that of hard 


spheres, the mean value of (l - cos ©) is unity, so that and 

are equal. This corresponds to an angular distribution function I q ( & ) 

independent of & . For very slow electrons this is a good approxi- 
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mation . In this case the collision frequency for momentum transfer 
may be defined as 



N Q v = N Q v 
m o 


( 68 ) 
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since the mean free path is equal to the reciprocal of the molecular 

number density times the cross-section, v is the mean velocity and 

l/p 

for a Maxwellian distribution is (8kT / ran) ' • 

The molecular number density is usually given in terms of 
the pressure normalized to 0°C or 


~] O 

N = 3.54 x 10 p (particles/cnr ) 


( 69 ) 


where v is in mm of mercury and is equal to 
- o 



Thus using equations 68 , 69 and 70> the cross-section 
for momentum transfer may be found from the collision frequency. 
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PART III. 

APPARATUS AND METHOD OF PROCEDURE 
A. The Plasma Cell 

In previous microwave transmission method experiments, 
a standard X-band waveguide has been tapered into a square section 
and a cylindrical glass tube with tapered ends has been placed in 
this square section to form the plasma cell. ’ The discharge 
takes place between electrodes placed external to the waveguide on 
each end of the glass cell. The positive column of a glow discharge 
is formed in the cylindrical glass tube within the waveguide. 

In the present experiment, the cell consists of a length 
of OFHC X-band waveguide with bakeable windows on each end. The 
entire cell is constructed of ultra-high vacuum materials suitable 
for bakeout to A00° centigrade. Holes are drilled in the top of the 
waveguide and anode wires extend down to these holes flush with the 
top of the waveguide. The waveguide acts as a hollow cathode and the 
negative glow of the discharge spreads down into the cell and out to 
several inches from the anode. The cell has two anodes approximately 
seven inches apart in order to provide a longer interaction path for 
the microwave signal. 

The advantages of the cell are twofold. First, the ques- 
tionable approximations and assumptions necessary with the enclosed 
glass cylinder mentioned previously are removed. The tapered glass 
cylinder causes some question as to the effective length of the micro- 
wave -plasma interaction and accounting for the presence of the glass 
involves approximations in obtaining the wavelength of the microwave 
signal in the cell. The rectangular cross-section of the waveguide 
cell provides a simple shape for calculation of the microwave -plasma 
interaction. The second advantage is the fact that the electron 
temperature in the hollow cathode negative glow does not rise as high 
as the electron temperature in a positive column. The electron tem- 
perature would then be expected to drop to room temperature faster 
in the negative glow, giving a longer period for obtaining data in an 
isothermal plasma. 
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The disadvantage of the cell is that the pressure range of 
operation is limited by the spread of the discharge into the wave- 
guide. In nitrogen, for example, discharges were obtained in the 
range 0.2 to about 7*5 millimeters of mercury. In helium, discharges 
were obtained from 0-5 up to approximately 30 or 40 millimeters of 
mercury. 

The pressure range might be increased by increasing the 
size of the anode holes in the top of the waveguide. This has the 
adverse effect of increasing the reflections of the microwave sig- 
nal, however. Some compromise must be made. At present, the hole 
size is 3/16" and the pressure range of operation is that mentioned 
in the previous paragraph. The voltage standing wave ratio of a 
single 3/16" hole in X-band waveguide is less than I.05. The VSWR 
of the complete cell is approximately 1.01. In order to obtain such 
a low VSWR, it was necessary to take swept frequency measurements 
on the cell and to pick an optimum frequency of operation. 

B. The Vacuum System 

A block diagram of the vacuum system is shown in Figure 4. 

All parts included in the "oven limits" are constructed of suitable 
ultra -high vacuum materials and are bakeable to at least 400° centi- 
grade. This section of the vacuum system is isolated from the remainder 
of the system by activated alumina traps which have a high affinity 

for the low vapor pressure organic oils used in the diffusion pump 
27 

and manometer. 

Measurements of pressure were obtained during roughing 
with a Veeco Type DVM1 thermocouple gauge. The system was baked at 
approximately 400° centigrade for approximately twenty-four hours 
while pumping with an oil diffusion pump (Consolidated Electrodynamic 
Corp. Type MCF-6l) to drive absorbed and adsorbed gases from the walls 
of the system. Following the baking, the metal ultra-high vacuum 
valves, numbered 1 and 4 in Figure 3, were closed, and the closed 
system was pumped to ultra-high vacuum pressures by an Ultek Series 
150 vac-ion pump. . The pressures were measured by a Veeco RG75 ion 
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gauge which had a low range of 10 mm of mercury, and since this 

range was achieved, the ultimate vacuum reached is not known, hut is 

-9 

at least 10 ' mm of mercury. The clean vacuum system leak rate is 

less than 10 ^ mm of mercury/minute. Since experiments were carried 

out at pressure in the vicinity of 1 mm of mercury, the impurity 

level due to vacuum system material contamination was less than 1 

8 

part in 10 . Gas samples used in all experiments were Linde Co- 
MSG (mass-spectrometer controlled) grade samples. These have a 
maximum impurity level of 10 parts per million for rare gases and 
100 parts per million for diatomic gases. These samples were leaked 
into the test cell through an ultra-high vacuum , valve . Pressures in 
the experimental range were measured with a low vapor pressure oil 
(Octoil-S) manometer isolated from the test section hy an activated 
alumina trap. 

C. The Microwave Bridge Circuit 

The microwave circuit used to obtain measurements of the 

phase shift and attenuation of a microwave signal passing through a 

decaying plasma is an extension of a basic microwave interferometer. 

In an interferometer circuit, a portion of the signal is coupled into 

a branch bypassing the plasma vessel and is then mixed with the signal 

passing through the plasma by use of a three port "T" junction. It is 

described adequately, for example, in S. C. Brown's book on Basic Data 
l8 

of Plasma Physics . The, interferometer pattern caused by a decaying 
plasma goes through maxima and minima corresponding to the phase shift 
of the microwave signal passing through the decaying plasma. The am- 
plitude of the fringe pattern is related to the attenuation of the 
same signal. 

The idea for the circuit used in these experiments came from 
an article by R. F. Whitmer^, who suggested that it might be possible 
to mix the interferometer signals in a hybrid T junction and obtain 
direct measurement of attenuation and phase shift from the two output 
arms of the hybrid T. 

That this is indeed the case may be seen by considering the 
properties of such a junction shown in Figure 5* Assume the electric 







fields of the TE 1Q mode are in the directions of the dotted arrows 
drawn at faces of the ports. The properties of the junction are that 
the electric field entering port 1 is split into equal fields leaving 
ports a and b with a phase reversal in port b. The signal entering 
port 2 is split equally into signals leaving ports a. and b with no 
phase reversals. There is no coupling between ports 1 and 2. Thus 
the instantaneous values of the electric field intensities may be 
given by 


e = 

a 


e 2 + S 1 



(71) 


If it is, assumed that is the electric field intensity 
of a reference signal and e^ is this same reference signal after being 
attenuated and phase shifted by a decaying plasma (or any other cause), 
the fields may be expressed mathematically as 


e 


1 





(72) 


where A is the attenuation of the signal entering port 2 in nepers 
and <t> is the shift in phase of this signal from the reference signal. 
If equations 72 are substituted in equations 71> "the electric fields at 
the output arms may be obtained. 
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The magnitudes of e & and e^ are of interest, since 
nitude is the quantity detected by a crystal detector mount, 
magnitudes a!re expressed as E & and E^, it may be found that 
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A family of constant attenuation curves on a plot of E^ 

versus E may be obtained, by eliminating (ft from equations 74. This 
81 

is easily accomplished by squaring both sides of equations 74 and 
adding. The result is 
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(75) 


which is the equation of a circle with radius 

H.--V ("l + e - 24 ) & 

Eliminating the attenuation. A, from equations 74 results 
in a family of quart! cs for constant phase lines given by 

E **■ - 2E 2 E 2 + E 4 = E 2 cos 2 ^ (E 2 + E 2 ) - E \os 2 d ( 76 ) 
a ab b 1 r ^ a b 7 1 r v 1 

These results are shown graphically in Figure 6 . Thus if 
E and E are measured, it is possible to obtain the attenuation. A, 

£L D 

and phase shift, ^ , from Figure 6. Note that the characteristics 
of the hybrid T blend well to afterglow problems where the attenuation 
is low. 

The microwave circuit is shown schematically in Figure 7* 

The circuit is initially balanced (E^ = 0) with the precision phase 
shifter in the lower arm and the attenuator in the upper arm. It is 
then possible to calibrate the bridge by employing the precision phase 
shifter and precision attenuator to plot a characteristic similar to 
Figure 6. In practice, the Figure 6 characteristic will be distorted 
somewhat by the nonlinear characteristics of the crystals and the non- 
ideal properties of the hybrid T detector. An actual calibration is 
shown in Figure 8. 
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Fig. 6. Hybrid T characteristics. 








Fig. 8. Microwave bridge calibration. 


One other consideration of the microwave circuit which must 
he taken into account is the heating of the electron gas by the alter- 
nating electric field of the microwave signal. Since the microwave 
signal is used as a diagnostic tool, its influence on the electron 
motion should he negligible. A criterion has been arrived at by V. 

E. Golant 81 by consideration of the change in mean electron energy due 
to the microwave field relative to the change in mean electron energy 
due to collisions with heavy particles. If the microwave field is to 


have a negligible effect, then the criterion may be stated mathematically 


as 
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( 77 ) 


where the term on the left represents the mean change in energy of an 

electron due to the microwave field and the right hand side represents 

the mean change in energy due to collisions with heavy particles. K g 

is the effective energy transfer coefficient per collision, W is the 

mean electron energy and is the mean energy of the heavy particles. 

If the real part of the conductivity is substituted in equation 77 and 

K is taken as 10 8 for a molecular gas^, at 300° K the criterion 
e 

reduces to 

E 2.1 x. 10" 8 (to 2 + 'V 2 m ) 1//2 volts/meter (78) 

The experimental conditions are that 'iS « <x> and a fre- 
quency of 9180 megacycles was used. If these conditions are used in 
equation 78, it is found that the peak electric field intensity of 
the microwave signal should be much less than 1210 volts per meter. 

In X-band waveguide operating at a frequency of 9180 megacycles, this 
peak field corresponds to a power of 157 milliwatts for the TE mode 
of operation. The maximum microwave power passing through the plasma 
cell in any of the experiments was 5 milliwatts, which is less than 
5/0 of this critical value. 

D. Gated Radiometer - Temperature Measurements 
The device employed for measurements of temperature in the 



afterglow is a modified version of a Dicke comparison type radiometer 
Basically, the comparison type radiometer compares radiation from two 
sources and, through a null detection system, indicates when they are 
equal. If one source of radiation is a Black -Body radiating at a known 
variaBle temperature and the second source is a plasma, the radiometer 
may Be employed to indicate when the two radiated powers are equal. As 
indicated in part II d, the radiation from a plasma is B(co,T) A , so it 
is then theoretically possible to calculate the temperature of the 
plasma if the absorption term A^ is known. It was also pointed out 
in part II d (equation 62), that if the plasma were located Between 
the standard and the radiometer, it would Be possible to determine the 
plasma temperature without knowledge of the absorption term. Since 
an afterglow eventually disappears and the radiometer would then re- 
ceive the radiation directly from the standard source, opportunity 
exists for switching the radiometer from the plasma to the standard 
By simply using a repetitive discharge of such a duration that the 
plasma decays completely in the first half of the discharge period. 

The details of this arrangement may Be derived from consideration of 
Figures 9 and 10 . Figure 9 is the microwave circuit consisting of 
(from the left) a termination; an argon discharge noise standard; a 
precision attenuator; the plasma cell; and the microwave receiver. 

The noise standard is essentially a Black -Body radiator at a known 

temperature T • The precision attenuator is at room temperature, 
s 

T q , and its attenuation, a , is accurately known. In the following 
discussion a is the attenuation ratio; i.e., a is equal to unity when 
there is no attenuation and to infinity when there is perfect ab- 
sorption. Thus, the absorption term, A , in equations 59 - 62 may Be 
expressed as 


A 


co 


= 1 



(T9) 


If the Black -Body radiation term B(oo,T) is expressed as its equivalent 
temperature, T, then By equation 60, the temperature equivalent of the 
microwave power at point T 1 in Figure 9 is 
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Pig. 9- Microwave circuit for temperature measurements. 



Fig. 10. 


Temperature measurements timing diagram. 
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Since all quantities on the right hand side of equation 80 are known 
and a is adjustable, then T-^ is a known variable temperature. 

If the equivalent temperature of the electron gas in the 
plasma cell is T g and the cell attenuation is a , the temperature 
equivalent of the microwave power "seen" by the receiver, T„, is 

X\ 

given by 
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(81) 


The variation in time of T_ is sketched in Figure 10. 

I\ 

The temperature and attenuation of the plasma cell rise rapidly during 
a short (5 microsecond) high voltage discharge (at T^ in Figure 10). 
The cell temperature, T , decays toward room temperature and by 
equation 8l, T^ is approximately equal to T g as long as the attenu- 
ation, a , remains appreciable. Then as the plasma disappears by 
recombination, etc., approaches unity and T^ rises to the known 

standard level T, • At the point where T decays through the Jknown 
_L xi 

temperature H } equation 8l becomes 

(1 - T e (1 - -5— > < 82 > 

e e 


Thus at this point T R is equal to T g and is also equal to 

the known value T^, irrespective of the value of the cell attenuation. 

The discharge repetition rate is 100 cycles per second, which 

allows sufficient time for the plasma to decay in the first half of the 

period between discharges. The IF of the microwave receiver is gated 

on at a 200 cps rate at the times designated as t in Figure 10. The 

s 

detected signal is the cross hatched pulses shown in Figure 10. The 
fundamental component of this pulse train is amplified in a high gain, 
narrow band, 200 cps amplifier and the result detected in a coherent 
phase detector. 

The time between the discharge t^ and the sample t is 

JJ s 
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adjustable; thus measurements are made simply by setting t at a known 

s 

time following the discharge and then adjusting the precision attenu- 
ator, a , until the microwave receiver indicates a null. T is then 
p' e 

equal to and is given by equation 80 . 

The gate width of the receiver is adjustable with a minimum 
of 1 microsecond. The temperature sensitivity is approximately 10 
degrees near room temperature , on the assumption that the plasma cell 
has an attenuation of at least 10 Db. during the instant of measure- 
ment. Meaningful measurements have been obtained with cell attenuations 
as low as 0.2 Db. 

E. Spectrograph! c and Light Intensity Measurements 

In addition to results obtained from microwave measurements 
of electron density and collision frequency, and measurements of the 
electron temperature, useful information may be obtained from the visible 
spectrum of the plasma. Both spectrographic measurements and total 
light intensity measurements were made. 

Spectrographic measurements were made in nitrogen with a 
3 l/4" by 4 l/4" Hilger quartz prism spectrograph. Originally it was 
intended to use the results of the spectrographic measurements to ob- 
tain information concerning the presence of the various molecular 
species of nitrogen, N, Ng, N^, and N^_. The intent was to separate 
the bands of the different species by means of a Bausch and Lomb mono- 
chromator (after band location had been established spectroscopically) 
and to read out the band intensities of the separate species as a 
function of time in the afterglow by means of a photomultiplier tube. 
Experimentally, this measurement is not too difficult and preliminary 
results were obtained on the bands of N^- This method was not pursued, 
however, since only bands attributed to N^ were obtained. No lines or 
bands attributed to N, N^ or N^ were found. This method might be 
employed in other gases, where more than one molecular species is 
present. 

Total light intensity measurements may be used to extend the 
range of electron densities obtained from the microwave method. In 
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the late afterglow, the radiation in the visible spectrum will be due 
to recombination of electrons and positive ions. Therefore the .inten- 
sity, I, is proportional to the product of the electron and ion den- 
sities. If only one type of positive ion is present, the assumption 
of charge neutrality requires the electron density to equal the ion 
density, thus 

I , a n e 2 ( 83 ) 

The light intensity measurements are more sensitive than the microwave 
method. Thus the electron density-time data obtained from microwave 
measurements may be extended, using a photomultiplier to convert the 
light intensity to a voltage signal and applying this to an oscillo- 
scope. The oscilloscope deflection will be proportional to the square 
of the electron density. The proportionality constant is obtainable 
from a calibration of the oscilloscope deflection at sufficiently high 
electron densities where the density may be obtained also by a micro- 
wave measurement. This method was successfully used to follow the 
electron density decay in the late afterglow. Because of the cali- 
bration problem and the fact that sufficient data could be obtained 
from the microwave interaction, this method was not pursued in the 
nitrogen experiments. 

F. Outline of Experimental Procedure 

The chronological order of experimental procedure and. the 
description of experimental apparatus not covered elsewhere are des- 
cribed in this section. The first consideration is to use as pure a 
gas as possible. Therefore, before each experiment was begun, a new 
one liter flask of Linde MSC grade gas was attached to the vacuum 
system. The vacuum system was then baked at approximately 400° centi- 
grade for at least 2k hours while pimping with a diffusion pump. 
Following the baking period, the closed system was pumped with a vac- 
ion pump to ultra-high vacuum pressures. During this pumping period 
the microwave hardware was assembled and the microwave bridge cali- 


52 



brated. If the vacuum system pressure rise was sufficiently slow in 
the absence of pumping to indicate no leakage and a "clean" system, 
the gas supply bottle was opened. The gas was then leaked into the 
vacuum system through a metal ultra-high vacuum valve. 

In the calibration of the microwave bridge, the deflection 
of an oscilloscope signal corresponding to the detected output. arm 
signals of the hybrid T detector is calibrated by the precision phase 
shifter and attenuator. A plot similar to Figure 8 is obtained where 
the -units of the vertical and horizontal axes are "centimeters of 
oscilloscope deflection" . 

The plasma is formed in the plasma cell by means of a hard 
tube pulser. The total -current is approximately 7 amperes and the 
duration is approximately 5 microseconds. Since the voltage drop in 
the cell is approximately 280 volts, approximately 10 millijoules of 
energy is imparted to the plasma in each breakdown pulse. The experi- 
ment is repeated every 10 milliseconds. Photographs are taken of the 
two hybrid T output arm signals and using the calibration curve, 
attenuation and phase shift of the microwave signal as a function of 
time in the afterglow is obtained. Electron density and collision 
frequency calculations were then performed on the IBM 650 computer 
using both equations 51 and 58 * In equations 51 the assumption of 
a uniform spatial distribution of electron density was used, and in 
equations 58 the spatial distribution was assumed to be sinusoidal, 
corresponding to a diffusion controlled electron density decay. The 
data corresponding to equations 58 was used with equation kk to deter- 
mine the recombination coefficient and ambipolar diffusion coefficient. 
The data corresponding to equations 51 was used in the early afterglow, 
where recombination is predominant, to obtain the cross-section for 
momentum transfer. 

The temperature measurements were carried out following the 
phase shift and attenuation measurements, since it was necessary to 
dismantle the microwave bridge and to set up the radiometer circuit 
shown in Figure 9* Care was taken to provide the same conditions as 
in the previous experiments. The pressure of the gas in the cell was 
reset to the values used in the bridge measurements. The hard tube 
pulser was also adjusted to give the same discharge conditions. 



PART IV. 

EXPERIMENTAL RESULTS 


A. Introduction 

In this section the results of measurements are presented. 
First, an amhi polar diffusion coefficient was obtained in helium. 

Since this gas has been studied in some detail in several laboratories, 
a check on the experimental method and apparatus is obtained. Detailed 
results of temperature and electron loss coefficient measurements in 
nitrogen are then presented. Finally, preliminary results of measure- 
ments of the cross-section for momentum transfer are presented. Dis- 
cussion of results is deferred to part V. 

B. Ambipolar Diffusion Coefficient in Helium 

In order to obtain a check on the experimental method and 

apparatus, measurements of the ambipolar diffusion coefficient in 

helium were made. As stated above, helium gas has been studied in 

some detail in several laboratories 1 ^' ^ . 

The initial results for the ambipolar diffusion coefficient 

times pressure, D p, were scattered in the range 700 - 900 cm 2 Torr/sec. 

£1 

A spectrograph! c study also showed the presence of an impurity taking 
an active part in the discharge. It might be pointed out that other 
investigators have also found it very difficult to obtain helium dis- 
charges free from impurity ions, owing to the high ionization poten- 
tial of helium relative to the ionization potential of common impurities. 

A cataphoresis pump^ 1 was added to the vacuum system in 
order to further purify the helium. The resultant value of D p was 
730 + 50 chi Torr/sec, which compares favorably with the value ob- 
tained by Kerr and Leffel^ (700 cm 2 Torr/sec) in 1959 and the value 
obtained by Oskam ^ (665 + 25 cm 2 Torr/sec) in 1957* In all of these 

cases, the helium molecular ion, He 0 + , is expected to be the predomi- 

2 48 

nant ion taking part in the ambipolar diffusion in the afterglow . 

C. Results in Nitrogen 
1. Temperature Measurements 

The results of the temperature decay at various pressures 
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are shown in Figures 11 - 17- In all instances the time corresponds 
to the center of the receiver gate. The receiver gate widths used 
in obtaining each point are indicated at the top of each graph. The 
vertical limits on the plotted points give an indication of the sensi- 
tivity of the radiometer measurements. The standard noise source used 
was a TD-18 Bendix tube which has 15.83 Db excess noise, corresponding 
to a temperature of 11,390 degrees Kelvin. 

Over the entire pressure range, the maximum time required 
for the electron temperature to relax to room temperature is 250 
microseconds. The data used to obtain electron loss coefficients was 
therefore taken only at times longer than 250 microseconds. 

In order to facilitate discussion in part V, a composite 
graph of the temperature decay for several pressures is plotted in 
Figure 18. In order to reduce confusion, the points which have a 
vertical range in Figures 11 - 17 have been plotted as a discrete 
point using the center temperature of the vertical range. Also, the 
abscissa is limited to the first 150 microseconds of the afterglow, 
since all curves are very close to each other following this time. 
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Pig. 11. Temperature decay at 7.1 Torr. 
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Fig. 12. Temperature decay at 4.6 Torr. 
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Fig. 13. Temperature decay at 2.95 Torr. 
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Fig. l4. Temperature decay at 1.95 Torr. 



Time - microseconds 

Fig. 15- Temperature decay at 1.2 Torr. 








2. Electron Loss Coefficients 

It might be advantageous at this point to summarize the 
steps taken to obtain the ambipolar diffusion coefficient and recombi- 
nation coefficient. Photographs were taken of the oscillograph traces 
of the hybrid T detector output signals. The photographic negatives 
were placed in an enlarger so the curves could be traced directly 

onto graph paper. Thus E and E, (Figures 7 and 8 ) were obtained as 

Sj d 

a function of time. At specific times the attenuation and phase shift 
of the microwave signal were obtained from the calibration curve (Fig- 
ure 8 ). Time increments of 10 microseconds were used in the early 
afterglow where the electron density is changing rapidly and a 50 
microsecond time interval was used in the late afterglow. The. values 
of attenuation, phase shift, and the corresponding time were then 
tabulated and punched onto IBM cards. The electron density and col- 
lision frequency were calculated by machine (IBM 650 ) employing equa- 
tions 58 and the time corresponding to each calculation was trans- 
ferred to the answer card. The computer thus furnished tabulated 
values of electron density, collision frequency and time. 

The incremental changes in electron density and time were 
then computed. The computed results were used to obtain values for 
the left hand side of equation 44 



and the left hand side of equation 44 was plotted as a function of 
■n • k and a. were obtained from the intercept and slope of the resul- 
tant straight line by the method of least square deviation. D was 

8 , 

then found from k by the relation 


k 


D it 2 ( 
a ' 



+ 




As a check on the results, k and cc were substituted into 
equation 42. 
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(42) 


n o (t) = 


C^ke 
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1 - dC-je 


-kt 


where 


n (o) 
o v 7 

C 1 ~~ k + cai o (o) 

C was obtained from one of the data points and the value 
of n from equation 42 was plotted on the same graph as the data points. 
The resultant graphs axe shown in Figures 19 - 2 6. In these figures 
the solid line is n versus t from equation 42. The points are the 
actual data points and the dashed line is the density decay assuming 
loss hy the fundamental diffusion mode. The dashed line is the asymp- 
tote of the solid curve. The difference between the dashed line and 
solid curve is the effect of the recombination in equation 42. In 
the early period immediately following the discharge, it is expected 
that the solid curve will deviate from the data points due to losses 
in higher order diffusion modes and to changes in the recombination 
coefficient and the ambipolar diffusion coefficient, the loss co- 
efficients in general being temperature dependent. 

The values of the constants used in equation 42 and the 
ambipolar diffusion coefficients obtained are shown on each of the 
Figures 19 - 26. The results of the ambipolar diffusion coefficient 

times pressure, D p, are plotted as a function of pressure in Figure 
8- 2 
27. The average value of D p at pressures below 4.6 Torr is 123 cm 

-1 a 2 -1 
Torr sec with a spread of + 24 cm Torr sec . The point at 7-1 

Torr is extremely high, which may be due to the fact that at this 

high pressure diffusion loss is small compared to the recombination 

loss. The method employed to obtain the diffusion coefficient is 

relatively insensitive at this high pressure. 

The results for the recombination coefficient are shown in 

Figure 28. It is noted that the recombination coefficient increases 

with pressure initially, and tends to level out at high pressures at 

-7 -i _i 

approximately 7 x 10 cm J sec . These results will be discussed 
further in part V. 
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Fig. 19- Electron density as a function of time at 

7.1 Torr. 
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Fig. 22. Electron density as a function of time at 
1.2 Torr. 
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Electron density as a function of time at 
0.75 Torr. 
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3- Cross-section for Momentum Transfer 

Preliminary results have been obtained for the cross-section 
for momentum transfer. The results shown in Figure 29 were taken .at 
two values of pressure, 7-1 and 4.6 Torr. The points obtained at 7*1 
Torr are widely scattered, and the results obtained at 4.6 Torr are 
but slightly better. 

The method used was to calculate the cross-section from 
the measured value of collision frequency as outlined in part II E. 

The scatter of points can be attributed to the lack of sufficient 
accuracy in the attenuation measurements which were used to calculate 
the collision frequencies. The primary concern was to obtain accurate 
electron densities in order to calculate electron loss process co- 
efficients. The electron density depends primarily on the phase shift 
and it is insensitive to changes in attenuation in the late afterglow. 
Therefore, no attempt was made to obtain attenuation more accurately 
than the values obtainable from the hybrid T characteristic (Figure 8). 
In equation 51? however, it is noticed that the collision frequency 
is directly proportional to the attenuation. In order to reduce the 
scatter of points, more accurate measurements of attenuation would 
have to be made. All values shown in Figure 29 are within a factor 
of two of results obtained by other methods. 
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Fig. 29. Cross-section for momentum transfer. 
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PART V. 


DISCUSSION OP RESULTS 

The results for the recombination coefficient differ markedly 
from previous measurements using the guided transmission technique by 

Bialecke and Dougal" 1 "^. In their measurements reported as 300° K, "the 

-7 

recombination coefficient was found to be approximately 9'0 x 10 

O 

cmr/sec independent of pressure in the range 0.2 to 2 Torr. However , 
in their work they allowed only 50 microseconds for the electron tem- 
perature to relax to room temperature following the discharge, and the 
data were taken for only approximately 200 microseconds following the 
discharge (see Figure 2, reference 13)* On comparing these time values 
with the results of temperature measurements (Figures 11 - 17 ), we see 
that Bialecke and Dougal took their data when the electron temperature 
was much higher than room temperature. They also found that the recom- 
bination coefficient decreased with increasing temperature. One would 
expect then that their measured value of recombination coefficient would 
be lower than the value at 300°K- That this is not true may be ex- 
plained as follows. They used a linear plot of l/n g as a function of 

time and obtained the recombination coefficient from the slope, neglecting 

32 

electron loss by diffusion entirely. Gray and Kerr have done exten- 
sive work examining the reliability of results for the recombination 
coefficient based on l/ n e vs. time plots. They have concluded that 
one cannot deduce a reliable recombination coefficient from a plot of 
l/n^ vs. time, even though it may be linear, unless the range of elec- 
tron densities involved covers at least an order of magnitude, and that 
in general the method yields too high a value of the recombination co- 
efficient. Since Bialecke and Dougal obtained linear plots of l/n g vs. 
time for only a factor of three or four in electron densities, the 
results should be viewed with extreme caution. One additional points- 
Bialecke and Dougal have attributed recombination to the dissociative 
recombination via Ng + according to the reaction 

+ _ * * 

N + e N + N 

2 

where the asterisk denotes excited atoms. Such a. reaction would result 
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in spectra of atomic nitrogen, tut as stated in the results 

section of this present work, only spectra attributed to the molecule 

N were found. In a spectrographic study by Bryan, Holt and Olden- 
^ 21 

berg , lines attributed to atomic nitrogen were never observed in the 
afterglow over the pressure range of 3 - 30 Torr. 

It is therefore important to point out that, although the 
microwave technique is a powerful diagnostic tool, other types of 
measurements are necessary to complete the description or to clarify 
the questionable assumptions of afterglow studies. 

Other measurements of the recombination coefficient in nit- 
rogen were the following. Biondi and Brown‘d made the first microwave 
determinations in 1949* using a cavity method, and obtained 14 x 10 ^ 

_ r 7 O 

to 18 x 10 cirr/sec. at pressures between 5 and 10 Torr. Their data 

28 

were not corrected for diffusion loss. Faire and Champion in 1959 

obtained a value of 4 x 10 cur 5 / sec. at 4 Torr, which is in good 

agreement with the present work. They used the cavity method and 

corrected their data for diffusion loss. A value of l4 x 10 ^ cm^/sec. 

was obtained before corrections for diffusion were made. This is 

precisely the value obtained by Biondi and Brown. Quite recently 

34 

Kasner, Rogers, and Biondi ^ used a cavity technique and, in addition, 
a Boyd type rf mass spectrometer for ion identification. They 
also used an inert buffer gas (helium or neon) to inhibit particle 
loss by diffusion. In this way they were able to measure recombination 
coefficients at low partial pressures of nitrogen, where ordinarily 
diffusion would be the principal electron loss mechanism. The value 
of recombination coefficient found for N p + at partial nitrogen pres- 
sures of approximately .005 Torr was 5-9 + 1 x 10 cnr/sec. a sig- 
nificant result of their work is the fact that at nitrogen pressures 

of 0.1 to 7*0 Torr the predominant afterglow ions were N + and N< + . 

^ 

It was only at nitrogen pressures less than .01 Torr that was found 
to be the only significant ion in the afterglow. This fact will be 
discussed further in attempting to describe the afterglow process. 
Possible recombination processes in nitrogen are: 
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(1) N 2 + + e" — ► Wg + h'V 

(2) N + + e“ N + ht/ 

(3) N + + e“ + Ng — N* + N 2 

(4) Ng + Kg + N 4 + 

N u + + e" — ^ Ng + Ng* 

(5) Ng + + e" + Ng Ng + Ng* 

(6) Ng + + e" — >■ N* + N* 


radiative 

radiative 

3-body colli sional 

2- step dissociative 

3- body collisional 
dissociative 


Reactions 2, 3 and 6 are eliminated because the spectrum 

of atomic nitrogen could not be found in the afterglow. In an article 

45 

by Massey in which extensive theoretical work on recombination is 
reviewed, it is shown that the radiative recombination coefficient is 
expected to be less than 10 cnr/sec. Therefore reaction 1 above 
cannot produce the results of the present experiment. In the same 

45 

paper the three-body recombination coefficient for air "might be 

“10 Q 

as large as 2 x 10 p cnr/sec." (p is the pressure in Torr). If the 

recombination coefficient for pure nitrogen is within an order of 

magnitude of this value, the results of the present experiment are 

three orders of magnitude higher than the expected value for three - 

body recombination. Therefore reaction 5 is not likely. 

9 

Bates has shown that dissociative recombination is likely 
to occur faster than other processes; recombination coefficients in 
the order of 10 cnr/sec. are made likely in this way. This is the 
reason that Bialecke and Dougal chose reaction 6 to explain their 
measured values. However, since only Ng spectra are found in the 
afterglow, it is more likely that reaction 4 is predominant. This 

theory is also substantiated by Biondi’s work with the mass spectro- 

+ ■+■ 
meter showing N^ ions outweighing the population of the Ng ion in 

the pressure range of these experiments. 
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Varney^ in 1953, while studying drift velocities in nitro- 
gen, obtained two distinctive regions in plots of drift velocity vs. 
E/p which he attributed to the presence of the ion at high. E/p 

values (low pressure) and N, + ion at low E/p values (high pressure). 

' 4 55 

This attribution was strengthened in 1957 when he studied the tem- 
perature dependence of the drift velocities. He found that higher 

■4* -f- 

temperatures cause onset of the transition from to at lower 
E/p values. This change is in accordance with the belief that as 
impacts between ions and gas . molecules grow more energetic, the N, 

■f - 

ions should dissociate and travel at the speed characteristic of 
ions. He has shown that the reversible reaction 


N 2 + + N 2 

V + N 2 W 2 + + 2N 2 


(84) 


favors the presence of at higher pressures and N 2 at lower pres- 
sures (no specific pressure values given). 

From these facts, the recombination process in this present 
work is believed to be predominantly 


N^ + + e" N 2 + N * (85) 

The rate of this reaction would be expected to depend on 
the pressure, since according to equation 84 the presence of N 2 mole- 
cules is required to form the ion. The question is not entirely 
answered, however, since the 3 -body reaction 


N g + + e - + N 2 N 2 * + (86) 

although considered unlikely by Massey, is not entirely ruled out. 
Spectra would again show only N 2 lines and the process is pressure 
dependent. The recombination process may be a mixture of the pro- 
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cesses described, toy equations 85 and 86 changing as a function of 

■f + 

pressure as the relative populations of and are varied. Per- 
haps an extensive mass spectrometer study of the afterglow ion currents 
would answer several questions, hut one must he careful in such measure- 
ments, since the binding energy of N, + is low and dissociation through 

^ 52 

collisions may take place easily in a mass spectrometer system. 

The ambipolar diffusion coefficient obtained toy the micro- 

wave transmission technique has not been previously reported. As 

2 

mentioned in the results section, D p was found to he 123 cm Torr 

-1 2 a -1 

sec with a spread of + 24 cm Torr sec . The ion mobility may 
be calculated from this value using equation 19> The mobility may 
toe compared with the results of other determinations if it is referred 
to the normal temperature and pressure values of 0° C. and 760 Torr. 

This may be done by the formula 

7.63 D p / 

a o cm/ sec 

^o T volt/ cm 

g ' 

2 "* 1 . o 

where Dp is in chl sec Torr and T is the gas temperature, 300° K 

a ° ® 2 - 1-1 

in this case. The resultant mobility is 2.85 cm sec volt . This 

53 

value compares favorably with 2»9 obtained by Tyndall and Pearce in 

28 

1935 by drift tube measurements. Faire and Champion deduced a mo- 

2 -1 -1 

bility at standard conditions of 2-8 cm sec volt from measure- 
ments of the ambipolar diffusion coefficient, employing the microwave 
cavity technique. 

These results, all in good agreement, raise a perplexing 

question, however. The mobility of N + in nitrogen would he expected 

7 ^ 46 

to show effects of charge transfer • Results of Mitchell and Ridler 

for mobilities of various gases in nitrogen in which charge transfer 

19 

would not be expected are shown in Figure 30- From theory one would 
expect that the mobility in such a situation would be nearly propor- 
tional to (l + Mg/M^) 1 ^. M is the mafes of the nitrogen molecule 
and M_^ is the mass of the ion. This trend is evident in Figure 30 . 
However, it is also noted that the value obtained in the present work 
falls on the curve. The effect of charge transfer has been observed 
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Mass of ion 


Fig. 30. 


Mobility in nitrogen of various ions as a function of mass 
at 1 atmosphere pressure. 



in rare gas ions diffusing in their parent gases in the following 
way. Mobilities of alkali ions in the inert gases were obtained. 

A plot of mobilities vs. mass of ions similar to Figure 30 was plotted, 
but it was found in all cases that the point corresponding to the 
inert gas ion would fall at a point 20 to 60 $ below the curve. This 
decrease in mobility was attributed to charge transfer. The theory 
for the charge exchange effect is not well understood, but some 
attempts- have been made to study this effect in helium . From re- 
sults of the present measurements it would appear that either charge 
exchange effects are not as important in the case of molecular ions, 
or that another process may speed up the diffusion, compensating for 
the decrease by charge exchange . 

The results of the temperature measurements are' interesting 
in that they reveal a source of energy in the early afterglow after 
the discharge source has been removed. Figure 18 shows the decay of 
temperature at several different pressures. It is noted that there 
is a rise in temperature reaching a maximum and then a final decay to 
room temperature * Not shown in Figure 18 is an initial rapid drop in 
temperature following the discharge that is difficult to follow be- 
cause it occurs in a time comparable to the minimum gate width of 
the microwave radiometer. The point to notice is that the maximum 
temperature value occurs later in the afterglow at lower pressures. 

Two possible explanations of high temperatures in the nitro- 

29 

gen afterglow have been suggested by Formato and Gilardini . First, 
"some nitrogen molecules are excited to metastable levels during the 
discharge ; metastable -metastable collisions in the afterglow may result 
in the ionization of one of them, so that a continuous supply of high 
energy electrons is present in the afterglow". Secondly, "the electron 
ion recombination coefficient may decrease with increasing electron 
energy; low energy electrons are then rapidly removed and the energy 
distribution is shifted upwards". 

The second suggestion seems to be ruled out, since one would 
expect that this would lead to a monotonically decreasing energy, and 
could not explain the increase in temperature. The first suggestion 
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of the presence of metastables appears to be more promising, but even 
this would require an increase in metastable -metastable collision rate 
after the discharge source is removed in order to obtain an increase 
in temperature. 

48 

Phelps and Brown have argued that in helium, the electron 
density actually increases owing to metastable -metastable collisions 
following high energy discharge pulses through a reaction 

, T m , TT m „ , „ + 

He + He — ►- He + He + e 

and the related reaction in which a molecular ion is formed 

TT id TT m . .. + - 

He + He — He^ + e 

It is possible, admittedly only through conjecture, that 
analagous reactions in nitrogen 

W 2 m + N 2 m _^ N 2 + W 2 + + e" 

N 2 m + N 2 m -^N 4 + + e' 

may cause high energy electrons to appear in the afterglow and to 
increase the electron temperature. The conjecture is not completely 
without basis, however. First let us consider a mechanism for in- 
creasing the population of metastables following the discharge. 

The equations describing the rate of change of the meta- 
stable density, M, and the electron density, n, are 

= -kM + ^(on 2 ) - PM 2 (86) 

= - an 2 + pM 2 (87) 

loss per metastable molecule due to destructive 


<)m 


and 


^ n 


k is the rate of 
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collisions with. neutral molecules in the volume of the gas and also 
may describe the loss represented by the fundamental diffusion mode. 

is the fraction of metastable molecules formed per recombination 
event. In equation 87 , it is assumed that the loss of electrons is 
primarily by recombination. PM 2 is the rate of formation of elec- 
trons by metastable -metastable collisions. Superelastic collisions 
between electrons and metastables may also provide energetic elec- 
trons in the afterglow. In a weakly ionized gas this effect would be 
small in comparison to the destruction of metastables by neutral par- 
ticles. 

Assume also that PM 2 is negligible with respect to the rate 

of destruction, kM, in equation 86. The recombination coefficient de- 

2 

creases rapidly with increasing temperature. Thus the formation 
term in equation 86 increases during the rapid initial decrease in 
electron energy in the first few microseconds of the afterglow. This 
would result in an increase in the metastable density. The metastable- 
metastable collision process would then have an increasing effect 
during the early afterglow. 

Secondly, the time retardation on the peak temperature with 
decreasing pressure may be explained as follows. At higher pressures 
the collision frequency for metastable -molecule collisions would be 
increased. Thus at lower pressures the destruction of metastable 
states through neutral collisions would proceed at a slower rate, 
providing for the longer effect noted in the temperature measurements. 

The problem in presenting a quantitative analysis of the 
metastable effect is that the rate of formation and destruction, and 
the initial density of metastable molecules are unknown. The tempera- 
ture dependence of the recombination coefficient is also not precisely 
known. In reference to equation 87, it is seen that all of these 
factors would have to be obtained in order to determine if the elec- 
tron density initially increases. 

As a further check on this process, measurements of elec- 
tron density in the very early afterglow are in progress in an attempt 
to ascertain if the electron density initially increases. 

The results of the preliminary measurements for the cross- 
814 - 



section for momentum transfer were given in Figure 29* Although the 
results are scattered owing to the inaccuracy of attenuation measure- 
ments, it is interesting to compare the results with other determi- 
nations. This is test done ty calculating the collision probability 

P from the cross-section where 
m 


P 

m 


N Q 


m 


The curves shown in the papers of Phelps, Fundingsland and 
Brown^ (microwave cavity), Anderson and Goldstein^ (microwave trans- 
mission), and Crompton and Sutton ^ (calculated from measurements of 
electron diffusion and mobility in dc electric fields) are reproduced 
in Figure 31 , along with the results from the present experiment. It 
is impossible to assign trends to the results of the present experi- 
ment, but, as mentioned in part IV 3 > the method, with increased 
accuracy of attenuation measurements, is expected to resolve this 
point, especially when combined with the direct measurement of elec- 
tron energy by the microwave radiometer. 
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PART VI. 

CONCLUSIONS AND SUMMARY 


A "brief summary of the experimental results is as follows. 
The room temperature recombination coefficient in the nitrogen after- 

m , r 7 O 

glow has been measured as 0.4 x 10 cnr/sec at 0.25 Torr rising to 
6.8 x 10 7 cm^/sec at 7*1 Torr. The predominant recombination mecha- 

4 “ 

nism is thought to be the dissociative recombination of the N^ mole- 
cular ion, 





although the three body reaction 



N, 




may be a significant process. The ion mobility at standard conditions 
of 760 Torr, 0°C, as calculated from the ambipolar diffusion coeffi- 
cient, is 2.85 cm sec volt , comparing favorably with previous 
drift tube measurements and microwave cavity measurements. The tem- 
perature measurements have indicated a source of energy in the early 
afterglow which has been attributed to the presence of metastable 
molecules. Preliminary results for the collision cross-section for 
momentum transfer are inconclusive, but the problems in the experi- 
mental procedure are resolved. 

The main conclusion reached is that the afterglow is of such 
a complexity that a combination of experimental procedures is required 
to separate and describe the various physical processes which are tailing 
place. The necessity for some of the questionable assumptions used by 
previous investigators has been eliminated by the - determination of the 
electron temperature by the gated radiometer and the simplification of 
the geometry for the microwave -plasma interaction. Still unexplained 
is the apparent absence of charge transfer effects on the mobility. 

Also the exact recombination mechanism is still unknown. It may be that 
an intensive study by mass-spectrometer and light spectroscopic tech- 
niques will answer some of these questions. 
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PART VIII. 


APPENDIX 

Approximations in the Solution of Equation 39> Page 19 

Equation 39; page 19; describes the time and spatial vari- 
ation of the electron density when loss occurs simultaneously by 
diffusion and recombination. 

—jf- = D a V 2n ~ ^ (39) 

(Al) 

Unfortunately, no exact solution is known and one is forced to make 

approximations in order to obtain the coefficients D and ct. In the 

a 

p 

late afterglow the usual assumptions made are that an is small and 

Ik 

that the fundamental diffusion mode predominates. Biondi and Brown 

28 

and Faire and Champion have used an approximating equation 

2 

= -kn - an (A2) 

where the constant k is the reciprocal of the time constant of the 
fundamental diffusion mode. In solving equation A 2, they assumed n 
to be spatially independent; hence the partial derivative becomes the 
total derivative in time. 

Redfield and Holt^ have used an approximation which is de- 
rived as follows: the spatial distribution of electrons is assumed 

to be that of the fundamental diffusion mode. In the rectangular 
waveguide cell this would be 


<)n 

c) t 


n = n (t) sin-: ^ - sin— ^ sin— — (A 3 ) 

o b b c ' 

If the product of the three sine terms is written as m(r) and equation 
A3 is substituted in Al, the result is 

— y-j: = -kn Q (t) - an o 2 (t) m(r) ( Ah ) 

Equation Ak is then integrated over the volume of the container- In 
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the waveguide cell the result is 


dn (t) 

— = -kn (t) - 

dt o v ' 



Chi 2 (t) 

o 


(A5) 


In the approximating equation (equation 40, page 20 ) used 
throughout this thesis, the spatial term m(r) in equation A4 was 
evaluated at the center of the cell where it has the value of unity. 
This results in the same form of equation as used by Biondi and Brown, 
and by Faire and Champion, so that comparisons of previous measure- 
ments in nitrogen could be made. Note that this would result in a 
value for the recombination coefficient which differs by a factor of 
8/jr or 0.258 from the method of approximation by Redfield and Holt. 
The difference factor, derived in a similar manner, for the two 
approximations in a cylindrical cavity is 0.275* 

The approximate solution used in this thesis (equation 42) 
is 


n (t) 
o 


C^ke 


-kt 


1 - QC^e 


-kt 


(A6) 


Some idea of the error involved may be found by substituting this 
solution into the original equation Al. If equation A1 is written in 
the form 


E = -If- - D a y 2 n + cm 2 (A7) 

E is zero for an exact solution. Now if equations A 3 and A 6 are 
substituted into equation A7, the time and spatial variation of the 
error expression, E, is 

r \ 2. 2 -2kt 

m(r)oC, k e 

g _ -L 

(1 - o^e ) 

where m(r) = sin- — - sin- ffi— sin— 

a b c 
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Since the error is a function of the spatial variables, the original 
assumption of the sinusoidal spatial variation is, of course, wrong.. 
Note that E is zero at the center of the cell and at the cell boun- 
daries, the latter being trivial since n is zero at the boundaries 
for all time. Since E decreases with time at all spatial positions, 
the distribution approaches one which is sinusoidal. The exponential 
rate is essentially twice the rate of the fundamental diffusion mode. 
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